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Introduction 


Internal  waves  are  a  ubiquitous  feature  of  the  deep  ocean,  occur¬ 
ring  on  spatial  scales  from  meters  to  kilometers  and  time  scales  from 
minutes  to  a  day.  A  unified  picture  of  oceanic  internal  waves  was  given 
by  the  kinematic  spectral  description  of  Garrett  and  Munk  [1972,  1979]. 
This  model  combined  linear  theory  and  experimental  data  to  yield  predic¬ 
tions  of  displacement  and  velocity  spectra  and  the  associated  coherence 
functions  for  both  moored  and  towed  sensors.  Since  they  are  best 
described  as  a  stochastic  phenomenon  with  a  wideband  frequency- 
wavenumber  spectrum,  statements  about  their  instantaneous  amplitudes  or 
velocities  have  little  meaning. 

By  contrast,  internal  waves  are  sometimes  observed  on  the  continen¬ 
tal  shelves  or  in  enclosed,  marginal  seas  as  quasi-deterministic  groups 
of  coherent  waves  with  well-defined  periods  and  wavelengths.  These  wave 
packets  generally  propagate  toward  land  from  a  localized  source  region, 
and  are  believed  to  originate  by  the  interaction  of  the  ocean  tides  with 
abrupt  changes  in  bathymetry.  The  waves  have  periods  of  minutes  to  an 
hour  and  wavelengths  ranging  from  50  m  to  a  few  km,  and  the  associated 
internal  displacements  range  from  5  to  as  much  as  100  m.  The  surface 
signatures  of  these  internal  waves  are  frequently  observed  in  either 
optical  or  radar  imagery  from  aircraft  of  earth  orbiting  satellites. 

Since  internal  waves  are  a  common  feature  of  the  oceans,  there  has 
been  some  interest  in  computing  their  electromagnetic  effects  in  order 
to  assess  their  role  in  the  oceanic  electromagnetic  environment.  Deep 
ocean  internal  waves  have  been  studied  by  Beal  and  Weaver  [1970]  and 
Podney  [1975]  using  simple,  two-layer  or  constant  buoyancy  frequency 
models.  Petersen  and  Poehls  [1982]  and  Chave  [1984]  produced  more  real¬ 
istic  spectral  descriptions  by  combining  the  Garrett-Munk  description  of 
the  internal  wave  field  with  an  appropriate  form  of  the  Maxwell  equa¬ 
tions  of  electromagnetism.  The  results  indicate  that  the  spectral  lev¬ 
els  are  fairly  high  in  the  ocean's  interior,  but  boundary  effects  reduce 
the  spectra  by  decades  at  the  seafloor  and  sea  surface.  Internal  wavo- 
induced  magnetic  fields  may  also  serve  as  a  significant  noise  source  in 
aeromagnetic  applications. 

This  report  extends  the  analysis  to  cover  the  electromagnetic 
effects  of  coastal  internal  wave  packets.  While  no  observations  of  the 
electromagnetic  fields  produced  by  them  have  been  reported,  the  large 
amplitude  of  shallow  water  internal  waves  is  expected  to  induce  substan¬ 
tial  electromagnetic  anomalies.  The  report  is  organized  into  five  sec¬ 
tions.  The  next  part  contains  a  review  of  oceanographic  observations  of 
internal  wave  packets  from  both  surface  ships  and  satellite  imagery. 
Section  3  outlines  the  hydrodynamic  theory  for  nonlinear  internal  waves. 
Section  4  contains  a  derivation  of  expressions  for  the  electromagnetic 
fields  produced  by  an  isolated  internal  soliton  both  at  the  seafloor  and 
above  the  sea  surface.  The  last  part  of  the  report  contains  a  discus¬ 
sion  of  the  theoretical  results,  including  a  detailed  look  at  the 
predicted  magnetic  fields  using  measured  purumoters  from  three  well- 
documented  occurrences  of  internal  wave  packets. 
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Introduction 


Internal  waves  are  a  ubiquitous  feature  of  the  deep  ocean,  occurring  on  spatial  scales  from 
meters  to  kilometers  and  time  scales  from  minutes  to  a  da\ .  A  unified  picture  of  oceanic  internal 
waves  was  given  by  the  kinematic  spectra!  description  of  Garrett  and  Munk  1972.  1979  .  This 
model  combined  linear  theory  and  experimental  data  to  \ield  predictions  of  displacement  and 
velocity  spectra  and  the  associated  coherence  functions  for  both  moored  and  towed  sensors.  Since 
they  are  best  described  as  a  stochastic  phenomenon  with  a  wideband  frequency-wavenumber  spec- 
truni.  statements  about  their  instantaneous  amplitudes  or  selociiies  have  little  meaning. 

By  contrast,  internal  waves  are  sometimes  observed  on  the  continental  shelves  or  in 
enclosed,  marginal  seas  as  quasi-deterministic  groups  of  coherent  waves  with  well-defined  periods 
and  wavelengths.  These  wave  packets  generally  propagate  toward  land  from  a  localized  source 
region,  and  are  believed  to  originate  by  the  interaction  of  the  ocean  tides  with  abrupt  changes  in 
bathymetry.  The  waves  have  periods  of  minutes  to  an  hour  and  wavelengths  ranging  from  50  m 
to  a  few  km,  and  the  associated  internal  displacements  range  from  5  to  as  much  as  100  m  The 
surface  signatures  of  these  internal  waves  arc  frequently  observed  in  either  optical  or  radar 
imagery  from  aircraft  of  earth  orbiting  satellites. 

Since  internal  waves  are  a  common  feature  of  the  oceans,  there  has  been  some  interest  in 
computing  their  electromagnetic  effects  in  order  to  assess  their  role  in  the  oceanic  electromagnetic 
environment.  Deep  ocean  internal  waves  have  been  studied  by  Beal  and  Weaver  1970  and  Pod- 
ney  1975  using  simple,  two-layer  or  constant  buoyancy  frequency  models.  Petersen  and  Poehls 
,1982  and  Chsve  1984  produced  more  realistic  spectral  descriptions  by  combining  the  Garrett- 
Munk  description  of  the  internal  wave  field  with  an  appropriate  form  of  the  Maxwell  equations  of 
electromagnetism.  The  results  indicate  that  the  spectral  levels  are  fairly  high  in  the  ocean’s  inte¬ 
rior.  but  boundary  effects  reduce  the  spectra  by  decades  at  the  seafloor  and  sea  surface.  Internal 
wave-induced  magnetic  fields  may  also  serve  as  a  significant  noise  source  in  aeromagnetic  applica¬ 


tions. 


ThiF  report  extends  the  analysis  to  cover  the  electromagnetic  effects  of  coastal  internal  wave 
packets.  While  no  observations  of  the  electromagnetic  fields  produced  by  them  have  been 
reported,  the  large  amplitude  of  shallow  water  internal  waves  is  expected  to  induce  substantial 
electromagnetic  anomalies  Tne  report  is  organized  into  five  sections.  The  next  part  contains  a 
review  of  oceanographic  observations  of  internal  wave  packets  from  both  surface  ships  and  satel¬ 
lite  imagery.  Section  3  outlines  the  hydrodynamic  theory  for  nonlinear  internal  waves.  Section  4 
contains  a  derivation  of  expressions  for  the  electromagnetic  fields  produced  by  an  isolated  interna! 
soliton  both  at  the  seafloor  and  above  the  sea  surface.  The  last  part  of  the  report  contains  a  dis¬ 
cussion  of  the  theoretical  results,  including  a  detailed  look  at  the  predicted  magnetic  fields  using 
measured  parameters  from  three  well-documented  occurrences  of  internal  wave  packets. 


Background 

Narrow-,  linear,  propagating  zones  of  disturbed  water  or  breaking  waves  in  an  otherwise  calm 
sea  have  been  reported  by  mariners  for  at  least  a  century,  especially  from  the  island  archipelagoes 
of  the  western  f’acific  and  Indian  Oceans.  These  are  undoubtedly  the  surface  manifestations  of 
very  large-scale  oceanic  internal  waves  The  first  scientific  evidence  connecting  these  surface 
effects  with  internal  waxes  is  due  to  Perry  and  Schimke  1965  .  who  made  repeated  bathythermo¬ 
graph  casts  through  several  lines  of  breakers  that  passed  the  research  vessel  Pioneer  in  the  Anda¬ 
man  Sea  north  of  Sumatra  Thex  were  able  to  infer  the  existence  of  internal  wave  packets  of 
80  m  amplitude  beneath  zones  of  surface  whitecaps  200~800  m  wide,  spaced  about  3200  m  apart, 
and  extending  from  horizon-to-horizon. 

Subsequently,  periodic  or  quasi-periodic  features  on  the  ocean’s  surface  were  detected  in  opt- 
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variationf  in  surface  water  velocity  that  sweep  together  surface  oils  in  sones  of  convergence, 
changing  the  optical  reflectivity  Ewing.  1950  .  The  second  mechanism  requires  focussing  of  capil¬ 
lary  wave  energy  in  zones  of  convergence  due  to  changes  in  the  surface  stress  field,  resulting  in 
increased  roughness  and  reduced  reflectivity  Gargett  and  Hughes.  1972  .  In  either  instance,  sub¬ 
surface  wave  structure  causes  surface  changes  that  may  be  sensed  remotely,  assuming  light  wind 
conditions  and  a  proper  solar  illumination  angle 

Similar  results  have  been  achieved  using  synthetic  aperture  radar  (SAR)  from  both  aircraft 
Elachi  and  .Apel.  1976  and  the  Seasat  satellite  Trask  and  Briscoe.  1983:  Apel  and  Gonzalez. 
1983:  Fu  and  Holt.  1984  .  Gonsistency  of  the  images  with  an  internal  wave  origin  have  been 
demonstrated  indirectly.  S.AR  offers  several  advantages  over  optical  imaging,  including  insensi¬ 
tivity  to  the  presence  of  clouds,  weather,  and  darkness  or  an  improper  angle  of  illumination. 

A  number  of  common  features  of  the  internal  wave  images  in  continental  shelf  regions  have 
come  out  of  these  studies.  First,  the  waves  occur  in  groups  or  packets  of  typical  width  3  to  5  km 
landward  of  the  continental  slope  The  packets  are  separated  b\  distances  of  order  15  or  80  km. 
corresponding  to  time  intervals  of  about  12  or  24  hours,  suggesting  a  connection  with  the  semidi¬ 
urnal  or  diurnal  tides.  Second,  the  wave  crests  are  oriented  parallel  to  the  local  bottom  topogra¬ 
phy.  and  may  be  refracted  significantiv  in  the  presence  of  abrupt  changes  in  bathymetry.  Third, 
the  wavelengths  within  a  given  packet  vary  from  200  to  4000  m.  depending  on  the  geographical 
area,  and  there  is  a  monotonic  decrease  in  wavelength  from  front  to  back  within  the  group. 
Fourth,  the  length  of  the  wave  crests  ranges  from  a  few  to  over  100  km.  and  the  crests  are  curved 
with  their  convex  sides  pointing  in  the  propagation  direction.  Finally,  as  the  packet  progresses 
onto  the  shelf  and  into  shallower  water,  there  is  evidence  for  an  increase  in  the  wavelengths,  sug¬ 
gesting  some  combination  of  nonlinear  and  dispersive  effects. 

The  behavior  of  large-scale  internal  wave  packets  have  been  observed  in-«itu  using  a  variety 
of  sensors,  including  moored  or  towed  thermistor  chains,  current  meters,  conductivity- 
temperature-depth  (CTD)  profilers,  expendable  bathythermographs  (XBT).  and  acoustic  scanners 
These  studies  may  be  classified  into  two  groups  according  to  the  water  depth  and  wave  ampli- 
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tude  The  shallow  water  and  continental  shelf  work  of  Halpern  1971  and  Haury  el  al.  1979  in 
Massachusetts  Hay.  Ingram  1978  in  the  St.  Lawrence  estuary.  Curtin  and  Mooers  1975  off  of 
the  Oregon  coast  and  Gargett  1970  near  \  ancouver  Island  is  confined  to  water  that  is  general!) 
less  than  ICHi  ni  deep  The  observations  typically  show  short  period  (5-10  minutes),  relatively 
low  amplitude  |5-Ki  ml  oscillations  of  the  thermocline  in  conjunction  with  surface  rips  or  linea- 
tions  of  the  type  seen  in  satellite  images.  The  packet  group  velocities  range  from  25  to  90  cm  s. 
with  the  large  spread  undoubtedly  influenced  by  advection  from  local  water  currents.  As  was 
inferred  from  satellite  sensors,  the  wavelengths  within  a  given  group  decrease  from  front  to  back, 
and  individual  elements  are  most  easily  distinguished  at  the  front  of  the  packet.  A  correlation  of 
the  occurrence  of  internal  waves  with  a  change  in  the  tide,  usually  from  ebb  to  flood,  is  notice¬ 
able.  and  the  spacing  of  wave  groups  is  also  consistent  with  a  tidal  generation  mechanism  For  a 
detailed  view  of  internal  wave  packet  evolution  in  shallow  water,  see  Haury  et  at.  1979  .  who 
used  high  frequency  sonar  to  view  their  development  in  time. 

The  deep  water  studies  of  Sandstrom  and  Elliot  1984  on  the  outer  continental  shelf  off 
Nova  Scotia.  Ziegenbein  19C9.  1970  near  the  Straits  of  Gibraltar.  Osborne  and  Burch  1980  in 
the  .Andaman  Sea.  and  Apel  el  al.  1985  in  the  Sulu  Sea  show  far  more  dramatic  internal  wave 
signatures.  Sandstrom  and  Elliot  1984  used  acoustic  and  CTD  methods  to  examine  internal 
waves  with  amplitudes  to  60  m  at  the  edge  of  the  Scotian  shelf.  Each  group  contained  only  a  few 
(2  4)  individual  elements,  and  were  seen  to  break  up  within  about  10  km  of  their  appearance 
near  the  shelf  edge,  a  behavior  at  variance  with  typical  satellite  observations.  They  also  sug¬ 
gested  that  mixing  by  these  tidally-produced  waves  was  a  primary  mechanism  to  supply  nutrients 
to  the  euphotic  zone  at  the  shelf-slope  break,  accounting  for  the  high  biological  productivity  of 
that  region.  Ziegenbein  1969.  1970  used  moored  thermistor  chains  in  water  800  to  1000  m  deep 
to  study  short  internal  wave  packets.  The  groups  of  internal  waves  were  associated  with  a  longer 
duration  internal  front  characterized  by  a  depression  of  the  thermocline  by  about  40  m.  Oscilla¬ 
tions  of  this  front  began  about  40  minutes  after  its  first  arrival,  with  amplitudes  of  as  much  as 
40  m  peak-to-peak.  The  group  velocities  of  the  wave  packets  ranged  from  1.6  to  2.2  m/s.  and  the 
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repeal  time  of  the  group?  indicated  a  lida!  a??oriation.  with  the  strongest  events  separated  b>  24 
hours 

Bn  far  the  best  documented  occurences  of  nonlinear  internal  vNave  packets  are  from  marginal 
seas  in  Southeast  Asia.  NN'here  a  combination  of  satellite  images  and  current  meter.  CTD.  XBT 
acoustic  depth  sounding,  and  surface  observations  have  been  applied.  Osborne  and  Burch  1980 
measured  packets  of  internal  waves  with  a  semidiurnal  recurrence  lime  that  was  markedly  modu¬ 
lated  by  the  tiNice  monthly  springtide  Amplitudes  to  60  m  were  seen  that  were  always  in  the 
form  of  depressions  from  a  normal  thermocline.  and  the  elements  of  a  wave  group  were  rank- 
ordered.  with  the  longest  wavelengths  at  the  front  The  ocean  stratification  could  be  roughly  fit 
by  a  two- layer  model  with  an  upper  region  of  230  m  thickness  underlain  by  an  860  m  thick  *one 
N'elocities  of  about  56  cm  s  in  the  direction  of  packet  propagation  occurred  in  the  top  layer,  while 
a  velocity  of  about  15  cm  s  with  the  opposite  direction  was  seen  below  this.  Surface  rips  were 
seen  above  the  descending  limb  of  the  internal  waves  cellular  flow  structure. 

.4pel  el  al.  1985  placed  three  current  meter  moorings  at  varying  distances  from  the  origin 
point  of  the  surface  lineations  seen  by  satellites  in  the  Sulu  Sea.  and  conducted  an  extensive  ship¬ 
board  profiling  program.  The  mooring  data  was  dominated  by  the  semidiurnal  tide,  but  large 
spikes  associated  with  internal  wave  packets  also  occurred  as  the  ebb  tide  turned  to  flood.  The 
amplitude  of  these  groups  was  heavily  modulated  by  the  14  day  spring  tide  cycle.  The  waves 
existed  as  depressions  of  the  thermocline  with  amplitudes  to  90  m  and  were  rank-ordered  and 
well-formed  at  the  leading  edge.  The  wave  trains  were  superimposed  on  a  broad  thermocline 
depression  of  about  50  m  amplitude  and  6  h  period.  Near-surface  velocities  of  up  to  1  m  s  were 
oriented  in  the  packet's  propagation  direction,  while  the  deeper  water  flowed  in  the  opposite 
direction,  and  the  turning  point  in  this  cellular  flow  was  at  about  200  m  depth.  Breaking  surface 
waves  were  seen  above  the  descending  limb  of  individual  packet  elements  with  amplitudes  of 
0.5-1  m  and  zone  widths  of  about  1  km.  The  mean  wavelength  within  each  group  was  2.6  km. 
and  coherent  wave  packets  were  seen  to  exist  for  up  to  2.5  d  from  creation  to  break-up  in  shallow 


water. 


A  generation  mechanism  for  interna!  wave  packets  by  tidal  flow  across  topographic  fealure^ 
has  been  investigated  theoretically  by  Lee  and  Beardsley  1974  and  using  laboratory  models  by 
Maxwc>rthy  1979  NN  hiie  they  differ  in  some  details,  the  overall  conclusions  of  these  studies  an 
similar  figure  1  outlines  the  basic  ideas.  As  the  ebb  tide  proceeds,  a  lee  wave  field  is  generated 
behind  the  topographic  ridge  by  the  locally  supercritical  flow.  The  magnitude  of  the  lee  wave 
field  depends  on  a  variety  of  factors,  including  the  shape  and  height  of  the  ridge,  the  water  velo¬ 
city.  and  the  stratificai ion  of  the  water  column.  As  the  ebb  tide  slackens  and  turns  to  flood,  the 
lee  wave  depression,  which  has  developed  a  substantial  phase  velocity  in  order  to  remain  station¬ 
ary  in  the  flt'w  ,  begins  tC'  move  upstream  against  the  tide.  As  the  front  continues  to  propagate 
away  from  the  ridge,  it  evolves  into  a  series  of  solitary  internal  waves  due  to  the  effects  of  non¬ 
linearity  and  dispersion.  The  si7e  and  number  of  these  waves  depends  in  a  compie.x  way  on  local 
properties  of  the  water  column  and  on  seafloor  topography 

Hydrodynamic  Theory 

The  observational  data  reviewed  earlier  indicate  that  internal  wave  packets  are  nonlinear; 
the  key  evidence  for  this  is  the  higher  velocity  of  the  longer  wavelength  components,  as  suggested 
by  the  presence  of  the  longest  wavelengths  at  the  leading  edge  of  any  given  group.  In  most 
instances,  the  individual  wavelengths  are  at  least  of  the  order  of  the  water  depth.  This  means 
that  the  appropriate  physics  is  that  for  nonlinear  waves  in  a  stratified  medium  under  the  long 
wave  approximation 

The  theory  for  nonlinear,  long,  internal  waves  is  reviewed  by  Osborne  and  Burch  1980  , 
.Apel  and  Gonzalez  1983  .  and  Fu  and  Holt  il984  .  For  a  two-dimensional  (i  ,r)  wave  field  (with 
z  positive  upwards  and  zero  at  the  sea  surface),  the  displacement  field  A  (*  ,z  ,1 )  may  be  written 

A  (t  .f  )  ^  *7(r  f  )0{z  )  (1) 

where  )  describes  the  horizontal  structure  and  ^(z)  is  a  dimensionless  vertical  mode  func¬ 

tion  which  is  normalized  so  that  its  maximum  value  is  unity.  For  waves  with  long  horizontal 
scales  L  and  small  (but  finite)  amplitudes  rf  in  water  of  depth  N  .  such  that  H  L  and  t)  //  are 
small  compared  to  unity,  perturbation  methods  e  g..  Benny.  1966;  Whitham.  1974  may  be  used 


10  get  relativelj  simple  equations  for  and  c  The  lowest  order  vertical  mode  equation  is  the 


familiar  linear  one  for  internal  waves 

S*©  -  ^  ~  ^  0  =  0  (2) 

ft* 

together  with 

0(0)  =  <?(-H  1  =  0  (3) 

where  -V(r)  is  the  Brunt-Vjiisala  or  buoyancy  profile  of  the  water  column  and  Co  is  the  eigen¬ 
value  of  (1).  which  corresponds  to  the  phase  velocity  of  long  linear  internal  waves.  This  equation 
IS  valid  for  frequencies  which  are  large  compared  to  the  local  inertia!  frequency .  but  small  relati\e 
to  the  buoyancy  frequency. 

The  zeroth  order  equation  for  the  horizontal  structure  function  is  the  Korteweg-deVries 
equation  Whitham.  1974 

d,T)  ~  -  at]d, Ti  ~  idit]  =  0  (4) 

The  environmental  parameters  a  and  *;  depend  on  the  buoyancy  profile  .V(r)  and  the  vertical 
eigenfunctions  o.  such  that 

0 

-Co  r  dz  d,.\'{:)0’{z) 

a  =  - ^ -  (5| 

r  dz 

-H 

0 

Co*  f  dz  <t>‘(z) 

-)  -  0  ■  ' -  (61 

2/  dz  y^(z]0^z) 

-H 

Apel  and  Gonzalez  [1983j  give  a  physical  interpretation  of  these  quantities.  The  nonlinear  param¬ 
eter  Q  is  a  measure  of  the  degree  of  amplitude  ’Stickiness"  or  the  degree  with  which  a  finite  ampli¬ 
tude  wave  travels  faster  than  its  nonlinear  counterpart.  This  is  caused  by  the  fact  that  a  finite 
amplitude  wave  perturbs  the  local  density  structure  significantly,  altering  the  local  dynamics. 
The  coefficient  q  is  a  measure  of  the  importance  of  physical  dispersion,  and  *  is  the 
wavenumber  at  which  a  quadratic  approximation  to  the  dispersion  relation  vanishes. 
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There  are  a  number  of  s[ead>-state  solution?  of  (4)  that  are  possible  when  the  nonlinear  and 
dispersive  terms  balance.  The  simplest  one  is  the  long-wavelength  internal  soliton  For  water 
that  is  deep  compared  to  the  vertical  scale  of  the  thermocline.  the  only  soliton  solution  that  is 
possible  IS  one  of  depression  given  by 


Vli  .t  }  =  -»?£,  sech*  ( — —) 

where  c  is  its  nonlinear  phase  speed  and  Z,  is  a  measure  of  its  width.  These  are  related  by 


r  =  c ,  1  - 


The  soliton  solution  (T)  consists  of  a  single,  bell-shaped  depression  of  the  thermocline.  and  is  a 
reasonable  appro.ximai ion  to  a  single  element  of  an  internal  wave  packet.  The  importance  of  the 
soliton  solution  lies  in  its  behavior,  including  1)  the  soliton  phase  speed  (8)  exceeds  the  linear 
phase  speed  c  ^  b>  an  amount  that  is  prof>orlional  to  its  amplitude  2)  since  larger  soliions 
travel  faster,  packets  of  them  will  evolve  into  rank-ordered  groups  with  the  largest  elements  at  the 
front.  3)  interacting  solitons  experience  at  most  a  phase  shift,  retaining  their  original  velocities 
and  shapes  after  collision,  and  4)  since  individual  solitons  evolve  with  independent  sizes  and  phase 
velocities,  the  distance  between  individual  soliions  in  a  packet  is  a  measure  of  the  distance  of  pro¬ 
pagation  from  the  source  These  properties  are  all  consistent  with  the  observed  characteristics  of 
internal  wave  packets.  Other  soliton  rules  are  reviewed  by  Osborne  and  Burch  1980  .  Note  that 
internal  wave  groups  may  be  regarded  as  superpositions  of  solitons  (7)  of  different  amplitudes,  and 
hence  different  velocities. 

A  second  steady-state  solution  of  (4)  is  reviewed  by  Apel  and  Gonzalez  1983  and  is  called 
the  cnoidal  solution  with 


viz  ,t)  =  f!m  -  (- 
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where  en„  is  one  of  the  Jacobi  elliptic  functions.  This  wave  also  propagates  without  change  of 
shape  and  consists  of  an  infinite  series  of  negative  troughs  of  amplitude  t][.  The  cnoidal  solution 
IS  a  nonlinear  generalization  of  the  cosine  and  reduces  to  that  form  in  the  linear  limit.  \\  hilt 
Apel  and  Gonzalez  1983  used  this  solution  to  model  observed  internal  wave  packets,  it  will  noi 
be  considered  further  in  this  report  due  to  its  analytic  complexity 

To  compute  the  environmental  parameters  (5)-(6).  it  is  necessary  to  have  a  good  estimate 
of  the  buoyancy  frequency  profile  .V(z  ).  This  may  not  be  available  in-iitu.  particularly  when  the 
data  consists  onh  of  satellite  imagery,  \^'hile  the  vertical  structure  equation  (2)  has  an  infinite 
number  of  eigensolutions  corresponding  to  the  different  vertical  modes,  observational  data  sug¬ 
gests  that  only  the  lowest  order  one.  having  a  single  maximum  at  mid-depths,  is  actually  present 
Ziegenbein.  1969.  1970;  Halpern.  1971;  Osborne  and  Burch.  1980.  This  means  that  numerical 
solution  of  (2)- (3)  for  the  first  mode  together  with  numerical  evaluation  of  (5)- (6)  can  allow 
realistic  approximations  to  be  made  with  the  internal  soliton  model,  commensurate  with  the  avail¬ 
able  knowledge  about  .V(r  )  This  procedure  was  used  by  Fu  and  Holt  1984  .  although  their 
buoyancy  profile  data  consisted  only  of  historical  averages  rather  than  synoptic  measurements. 

Instead  of  using  this  relativelj  complex  procedure,  and  in  view  of  the  many  other  approxi¬ 
mations  that  go  into  the  internal  soliton  model  (7).  the  simple  two-layer  model  of  Osborne  and 
Burch  1980  will  be  adopted  This  will  allow  analytical  expressions  for  the  electromagnetic  fields 
to  be  derived,  consistent  with  a  philosophy  of  simplicity  in  which  order-of-magnitude  estimates  for 
the  internal  wave  magnetic  fields  are  sought.  In  a  two-layer  ocean  model,  only  the  lowest  order 
vertical  mode  is  possible,  and  the  vertical  structure  function  in  (2)  is  approximated  by  straight 
line  segments.  For  an  upper  layer  of  thickness  h,  and  density  p,  overlying  a  lower  layer  of  thick¬ 
ness  and  density  pj.  where  h,  <  hj  so  that  (7)  is  valid,  the  environmental  parameters  (5)-(6) 
and  the  linear  phase  velocity  Co  are  given  by  [Osborne  and  Burch,  1980j 
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where  p  ~  ^1  ~  p;  and  Ap  =  pj-pj  The  horizontal  water  velocities  in  the  two  layers  are 

^  oVo 


(12) 


(IS) 


(H) 


a  ,(i  t  1 


A, 


sech'  $ 


,  .  <■  o’)  c  .  -.  , 

U;(r  .<  I  - 7 —  sech*  4 
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^  - 


I  -et 


(15) 

(16) 

(17) 


Note  that  the  horizontal  velocity  components  are  depth-independent  within  a  given  layer  but 
discontinuous  at  the  thermocline  z  -  -A,,  and  that  they  have  opposite  directions  in  the  different 
layers  The  vertical  velocity  components  follow  from  the  continuity  equation  and  the  rigid  lid 
approximation 


u  ,(i  .i  .<  )  =  -(d,  u,(i  .1  ))z 


2go’)oz 

f,,L 


sech*^  tftnh  ( 


U’j(t  ,r  .«  )  =  ~(d,  U2(»  0)  *W) 

2co’)o(^  ) 


A,Z, 


sech*^  tanh  ( 


(18) 


(19) 


Note  that  u  ,  =  u  j  at  the  thermocline  2  =  -A],  and  that  the  vertical  velocity  depends  linearly 
on  the  depth  with  a  maximum  value  at  the  interface.  The  water  motion  is  prograde  elliptical, 
with  the  upper  layer  water  moving  in  the  direction  of  wave  propagation,  descending  at  the  lead¬ 


ing  edge  of  the  soliton.  and  rising  at  its  trailing  edge.  A  small  surface  soliton  accompanies  the 
internal  one.  and  is  neglected  in  this  treatment.  Figure  2  illustrates  the  geometry  and  soliton 
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Y>ehavior  of  ihe  model 

Elertromaguptir  Theory 

The  framework  for  the  theory  of  electromagnetic  induction  used  in  this  report  is  similar  to 
that  in  Chave  I9h4  The  Maxwell  equations  in  the  quasistatic  limit,  in  which  the  magnetic 
effects  of  the  displacement  current  are  neglected,  govern  the  electromagnetic  fields 


B 

0 

(20) 

ti ,  B 

0 

(21) 

p  (T  £ 

(22) 

where  i  l^  f.e  elecinc  field  B  is  the  magnetic  induction,  ft  is  the  magnetic  permeability  of  free 
space  o  IS  the  electrical  conductivity,  and  is  a  specified,  impressed  source  electric  current  den¬ 
sity.  For  induction  processes  that  are  driven  by  moving  seawater,  the  latter  is  given  by 

r=ao{r>?)  (23) 

where  ft  c  the  conductivity  of  seawater.  7  is  the  water  velocity,  and  F  is  the  geomagnetic 
induction,  assumed  sourceless,  curl-free,  and  time  invariant.  Replacement  of  B  with  F  in  (23) 
linearizes  the  set  (20) -(23)  by  eliminating  the  effect  of  feedback  by  the  induced  magnetic  field  on 
the  velocity  field. 

Solutions  of  (20)-(23)  will  be  sought  using  potentials  in  a  Cartesian  coordinate  system.  By 
virtue  of  (20).  the  magnetic  induction  may  be  written 

B  =  V  X  (n/)  -  VaS,*  -  (24) 

where  II  and  '1'  are  scalar  functions  which  represent  the  toroidal  and  poloidal  magnetic  (TM  and 
PM)  parts  of  the  magnetic  field  and  V  *  is  the  horizontal  derivative  operator  d,£  -t-  The 

source  current  in  (22)- (23)  may  also  be  expressed  in  terms  of  scalar  functions 

=  =/  -  V*T  ^  V  MT/)  (2.^) 

where  T  and  T  satisfy 


V/?T  =  <to(Va  ^v,,  F,£  -I-  <ro(V  A  z'j  PA 


(26) 


-  c?  o(  .  h  '^h)^z  ~  ^  ^  ) * X 


E  -  p  c('«  • 


spatial  variations  of  F  have  been  neglected  in  (26)-(27).  consistent  with  the  study  of  small-scale 
oceanic  phenomena 

If  a  IS  assumed  to  be  a  function  of  depth,  uncoupled  partial  differential  equations  for  11  and 
flt  may  be  obtained  by  writing  the  electric  field  in  a  form  like  (25 1.  substituting  into  (21)-(22), 
writing  out  the  Cartesian  components,  and  applying  complex  variable  theory  to  the  result,  yield- 


?n  -  a'd  .{b  ,n  er  ]  ~  fiob,n  =  -  HZ.  -  //p9,(T  <t  \ 


-  /iffd.'I'  =  -pT 


£  =  -  -  *  n  -  E  .  p 


,  dj  n  ■>  -  T  p  -  "  >  (d,  tit  f  ’ 


Independence  of  the  TM  and  PM  modes  is  assured  because  the  usual  electromagnetic  boundar> 
conditions  at  horizontal  interfaces  are  decoupled. 

It  is  especialK  convenient  to  solve  (29)-(30)  by  expressing  the  lime  variable  (e""*'  conven¬ 
tion)  and  the  horizontal  coordinates  as  Fourier  transforms  using  the  standard  pair 

ff  dzdy  e  ■'“-'»)/  (z,») 

oc 

x 

=  ff  (32) 

(2tr  ) 

Green  functions  for  (29)- (30)  may  then  be  found  for  a  constant  depth  ocean  of  uniform  conduc¬ 


tivity.  where 


n(t,/,2)  =  -H  f  dz  'j,(z,2  ')5(r  ')  ~  fi  f  dz  'd,  g,(z,z  ')f(2  ') 


^(k,l,z)  =  f  dz'g*(r.2')f(z') 


(34) 
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and 
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-  t  /I  ;  €  C 
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^-^0  p/*M^  pPAf  •"* 'J  p  PM^  r-;' 

e  ^r\  '  t  e  ~  a>  c  —  a>i  a/  ^  c 

j, (,.,-)  - - ^ -  - (36) 


23o(i-RrRre 


The  induction  parameter  is 


3  =  \'k^~l'-i  u  fta  (r  ) 

with  given  when  <t  (c  )  -  (t  (,  while  the  modal  reflection  coefficients  are 
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ipK 
r.w  _  °  '■ 
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-  1 


(37) 


(38) 


Rl 


PM  _  ^0 


;9oA-l 


^oA  - 1 


(39) 


(40) 


Expressions  for  the  modal  response  functions  K  and  A  may  be  found  in  Chave  and  Cox  1982’  for 
a  layered  medium.  For  the  purposes  of  this  report,  only  a  halfspace  of  uniform  conductivity  o  , 
will  be  considered,  reducing  (38)-(39)  to 


TM  _  ^0^1 -^1^0 


Rl 


Rr 


3o^  1  ■*  0 


|go-/g. 


(41) 

(42) 


The  Internal  Soliton  Electromagnetic  Fields  -  Theory 

The  electromagnetic  fields  may  be  found  by  Fourier  transforming  (15)-(16),  (18)-(19),  and 


(26)-(28)  to  get  expressions  for  the  source  current  scalars,  evaluating  (S8)-(84)  to  get  the  elec¬ 
tromagnetic  scalar  functions  Tl  and  ft,  taking  the  inverse  Fourier  transform,  and  using  (24)  and 
(81).  Using  the  Fourier  transform  identity 


■ '\V 


f  dk  (‘‘'sech-’f  =  7ri'*A-csch(-^)t'*'' 

2 


^  1 ,  <*<> 


T,(r  .  /  )  ^  a  ^Fy  a  iL^kz  csch( — ^ — )« 


Ti(z.tl=  -racFya,L^Hz~B)cschl-j^)e‘^‘'  (45) 

f,(--.0  =  T  o  oL=a,csch(-^)e^'  IF.kz  -iF,}  (46) 

T;i--  .n=  7r«Tt,L-a,csch(-^)€‘‘''(^',  *■(-’ )  (4"l 

E,(c  .O"  TT^Tp/',  a,<rL*csch(-^~)«  (48) 

E2(:-<)  =  -T  a  oFy  a,kL^csch{-^'^)e‘^‘  (49) 

where  O]  =  Cgt/g  fi^  and  Oj- 

Rewriting  the  TM  mode  expression  (33)  by  integrating  the  second  term  by  parts,  it  is  easy 
to  show  that  the  induced  part  of  the  TM  mode  electric  field  and  the  TM  mode  magnetic  field 
vanish  if  E  — d,  T  -  0  and  T(  //  )  =  0.  Using  (44)- (45)  and  (48)  -  (49),  it  is  clear  that  these  con¬ 
ditions  hold,  and  the  TM  mode  need  not  be  considered  further.  From  (31).  there  is  a  remaining 
TM  mode  electric  component  given  by 


(44) 

(45) 

V  j 

.  ■-  -t 

(46) 

vv1 

■::y 

(47) 

(48) 

(49) 

t  "J 

=  (50) 

This  is  just  the  quasistatic  electric  field  required  to  ensure  that  no  electric  current  is  associated 


with  this  mode,  and  that  no  current  flows  in  the  i  -  z  plane. 

The  remaining,  PM  mode,  electromagnetic  problem  is  most  easily  solved  by  writing 
(46)-(47)  as  the  sum  of  parts  due  to  F,  and  F,  induction.  Taking  the  Fourier  transform  in  time, 
this  gives 

f,‘(r)  =  2ir^aoF,L\ikz  uch  {^^)6(w-kc  )  (51) 

Tj*(^  )  =  ~2»r  VqF,  L^Ojk  {z  +H)  csch  }6  {u' -  ke  )  (52) 


Ti  (i  )  =  Ht  ‘c  oF:  L'a  I  csch  ( — - — (-1-  ke 


Tf  o'l '  I  =  ~i  2t.‘ a  (jF.  L‘a'2  csch  (— ^ —  )^(u-  - kc 


where  <5(r  )  is  the  Dirac  delta  function.  Owing  to  the  linearity  of  the  MaAwell  equations,  the  two 
pans  of  the  PM  mode  scalar  "It  due  to  (51)-(52)  and  (53)-(54)  may  be  computed  separately  and 
then  combined 

The  interest  in  this  report  is  confined  to  the  magnetic  field  at  or  above  the  sea  surface 


0)  and  at  the  seafloor  (z  =  ~H)  From  (24).  there  are  four  Green  functions  that  are 
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Combining  (51)  (52)  with  (55)-(56),  evaluating  (84),  and  taking  the  inverse  Fourier 
transforms  of  the  i  and  f  variables  yields 


dk  e-*^« 


it„ch(-^)(f?/**^l)(o,.4,-a,>l,) 


),**(o)  =  — / dk 


*c.ch(-^)(f?/"-l)(a,/i,-Oj/l,) 


where  d,  =  s/k^-ike  no ,  in  (39) -(40)  and 


Wi-ih  "’h 


>»2=  [((/?0*2-l)-«i’"(^0*2^1)*  ”"*’)*‘'’"*'  -  («i“" 
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Applying  (24)  give^  the  i  -  and  i -components  of  the  magnetic  fields  at  the  sea  surface  To  get 
the  magnetic  fields  above  the  sea  surface,  note  that  in  air  and  in  the  absence  of  sources,  tj*  obeys 
a  one-dtmensional  homogeneous  Helmholtz  equation,  vanishes  as  r  —  oc.  and  is  specified  on 


2  =  0.  SO  that 

'!>(?)=  ♦(0)e"  ‘  '  (61) 

The  airborne  magnetic  fields  follow  by  taking  the  appropriate  derivatives  of  (59)  after  multiplying 
the  kernel  bj  c"  *  “.  yielding 
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e--it»csch(-^)(ff/^-l)(a,A,-QjAj) 
- = - j5-r -  (63) 

where  a  =  k  .  Note  that  (62)  reduces  to  the  r -derivative  of  (60)  when  r  -•  0. 

A  similar  procedure  gives  the  induction  part  of  the  magnetic  fields  at  and  above  the  sea 
surface 
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where 


The  magnetic  fields  at  the  seafloor  are  obtained  in  an  identical  manner.  For  induction  h\ 


the  Fj  geomagnetic  field,  they  are 
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and  for  induction  by  the  F,  geomagnetic  field,  they  are 
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where 
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D2  =  (l-e''’’*’)(l-/f/**«^^»''e'^'*‘)/^(, 

The  field  component  expressions  (62)~{63)  and  {66)“(7])  are  easily  evaluated  using  a  fast 
Fourier  transform  algorithm.  Care  must  be  taken  to  ensure  that  the  result  is  properly  normalized, 
so  that 

^Jdkc  f{k)^  /  Ik,  )e  /=  0,  ....  .V  -  1  (72) 

-oc  ^  ,=0 

Note  also  that  the  magnetic  fields  display  even  or  odd  symmetry  about  (  =  0.  and  can  also  be 
written  as  cosine  or  sine  transforms. 


The  Internal  Soliton  Electromagnetic  Fields  ~  Examples 

In  order  to  assess  the  size  and  spatial  scales  of  the  internal  soliton  magnetic  fields  both 
above  the  sea  surface  and  at  the  seafloor,  it  is  necessary  to  measure  or  estimate  the  hydrodynamic 
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quantities  in  ((i2)~(63|  and  (66)~'(T1)  Due  to  the  wide  variety  of  settings  in  which  these  iniernal 
wave  packets  are  observed,  and  the  concomitant  range  of  the  relevant  parameters,  three  examples 
were  selected  that  span  the  range  of  soliton  scales  that  seem  to  occur  in  nature  The  first  of  these 
is  from  the  Andaman  Sea  experiment  of  Osborne  and  Burch  198(1  .  in  which  in  situ  measurements 
of  the  water  velocity,  stratification,  and  soliton  length  scales  are  available.  The  even  larger  inter¬ 
nal  wave  packets  seen  in  the  Sulu  Sea  were  excluded  because  the  simple  shallow-water  soliton 
model  is  not  strictly  applicable  Apel  et  al..  1985  .  The  second  example  is  taken  from  the  satellite 
observations  of  Fu  and  Holt  1984  in  the  Gulf  of  California  While  no  in  tiiu  measurements  are 
available  at  this  site,  the  study  is  thorough,  the  water  depth  is  typical  of  the  outer  continental 
shelf,  and  it  is  located  at  mid-latitudes.  The  final  example  is  drawn  from  the  Massachusetts  Bay 
experiment  of  Halpern  1971  .  and  is  believed  to  be  representative  of  shallow  water  and  small  soli¬ 
ton  scales.  Table  1  lists  the  water  depths  and  hj  for  the  two-layer  soliton  model,  the  soliton 
amplitude  tjp.  the  soliton  half-width  L.  the  layer  water  velocities  a,  and  o  j.  and  the  nonlinear 
phase  speeds  that  have  been  either  mea.sured  or  estimated  for  these  three  cases.  Except  for  the 
Andaman  Sea  data,  these  values  are  quite  crude,  and  probably  not  accurate  to  better  than  a  fac¬ 
tor  of  2  to  4 

The  internal  soliton  magnetic  fields  also  depend  on  the  geomagnetic  field  components  F, 
and  F,  .  A  simple,  axially-geocentric  dipole  model  of  the  geomagnetic  field  was  employed 

Ff,  ~  Ffl  cos  A 

Fz  =  -2FosinA 

where  Fq  =  3.05  >  10‘  nT.  A  is  the  geographic  latitude  given  in  Table  1.  and  Fjs-  and  Fz  are  the 
north  and  vertical  geomagnetic  components.  The  errors  associated  with  neglecting  the  off-center 
nature  of  the  real  geomagnetic  field  and  ignoring  its  east-west  component  are  small  at  low-  to 
mid-latitudes.  Note  from  (62)~(63)  and  (68)~(69)  that  the  F,  field  is  the  component  of  the 
geomagnetic  field  in  the  direction  of  soliton  propagation.  This  means  that  a  soliton  will  generate 
the  largest  possible  magnetic  field  when  traveling  north,  while  one  that  is  moving  to  the  east  will 
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interact  only  with  the  vertical  geomagnetic  field,  yielding  a  smaller  anomaly  Note  also  that  the 
components  (66)~(67)  and  (70)“(71)  are  invariant  under  changes  in  the  direction  of  soliton 
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propagation 

Figures  3“6  show  the  four  parts  (62)  (631  and  (66)  (67)  of  the  magnetic  field  for  the  Ana- 
man  Sea  data  tn  Table  1  at  and  above  the  sea  surface  for  a  fixed  time  t  .  Since  (62)  and  (67)  can 
be  written  as  sine  and  (63)  and  (66)  as  cosine  transforms  by  taking  advantage  of  the  symmetries 
in  their  respectrive  kernels,  an  antisymmetric  form  is  expected  for  the  first  pair  of  fields  and  a 
symmetric  one  for  the  latter  pair.  This  reflects  the  nature  of  the  underlying  driving  electric 
currents,  and  can  be  understood  by  examining  their  spatial  distribution.  Transforming  (51)~(54) 
back  to  the  space  and  time  domains  gives 

Ti*(i'.2.t)  =  OoFjOjZ  sech^f  (74) 

T;‘(t.2.<)  =  -<7  oF,  OjIj )  sech'f  (75) 

T,'(i  .r  )  =  a  qFj  a  tanh  ^  (76) 

Tj'(i ,/ )  =  <To/'^  tanh  ^  (77) 

Applying  (25)  gives  the  driving  current  density,  which  has  a  singlet  component 

(78) 

The  two  layers  denoted  b)  the  subscript  may  be  considered  separately,  and  have  a  similar  but 
oppositely-directed  effect.  In  the  top  layer.  (74)  gives  an  antisymmetric  or  odd  source  function 
about  {  =  0  that  may  be  thought  of  as  due  to  a  pair  of  opposing  line  currents  on  either  side  of 
f  =  0.  The  horizontal  magnetic  field  from  such  a  current  pair  will  also  be  antisymmetric  and 
display  relative  peaks  above  each  line  source.  This  is  seen  in  Figure  3.  The  vertical  magnetic 
field  will  be  symmetric  and  have  its  largest  magnitude  between  the  two  line  sources  where  the 
individual  contributions  add  constructively.  This  is  seen  in  Figure  4.  The  second  part  of  the 
upper  layer  source  current  given  by  (76)  has  a  symmetric  or  even  functional  dependence  that  may 
be  thought  of  as  due  to  a  single  line  current  at  ^  =  0.  The  horizontal  magnetic  field  from  such  a 
current  is  symmetric,  while  the  vertical  component  is  antisymmetric,  accounting  for  the  behavior 
in  Figures  5  and  6.  The  contributions  from  the  lower  layer  are  similar  but  of  opposite  sign.  For 
an  observer  located  above  the  sea  surface,  the  upper  layer  is  dominant  due  to  its  proximity,  and 
the  oppositely-directed  anomaly  from  the  lower  layer  is  broader  and  of  lower  amplitude.  This  is 
seen  clearly  in  the  wings  that  are  apparent  in  Figure  4.  An  observer  at  the  seafloor  would  be 


more  sirongl>  afferleri  b>  th<-  I  'wer  la>er.  Nuie  that  the  actual  currents  in  (T4)~(771  with  (781 
are  extended  in  the  i -direction,  and  also  var>  in  the  z -direction  in  the  case  of  (74)“(7n). 

The  sohton  model  is  a  propagating  one.  as  seen  in  the  argument  in  the  complex 
exponential  in  (62)~  (fcS)  and  (66)  “(71).  For  an  airborne  observer,  the  sensor  velocitj  t  is  much 
larger  than  the  sohton  phase  vek>city  c  .  and  time  may  be  regarded  as  fixed  while  the  sensor 
moves  over  the  sohton.  For  a  seafloor  observer,  the  position  i  is  fixed,  and  the  time  progression 
of  the  soliton  will  be  seen 

Figures  7  and  8  show  the  horizontal  and  vertical  magnetic  fields  for  the  .\ndaman  Sea  soh¬ 
ton  propagating  to  the  north,  and  are  the  sums  of  Figures  S  and  5  or  4  and  6.  respectively.  The 
purely  symmetric  anomalies  of  Figures  3~6  have  now  vanished.  A  sohton  moving  to  the  easi 
would  produce  only  the  fields  of  Figures  5  and  6.  since  induction  by  the  north  geomagnetic  field 
does  not  occur  The  peak-to-peak  amplitudes  of  the  induced  magnetic  fields  range  from  0.5~ 
3  nT.  and  are  much  larger  in  Figures  7  and  8  because  the  geomagnetic  field  is  nearly  horisonial  at 
very  low  latitudes  (/",«  •.  Fi,  }.  The  anomaly  wavelength.*  reflect  the  horizontal  soliton  scale  L. 
in  this  case  about  1  kni.  The  magnetic  fields  are  largest  at  the  sea  surface,  and  attenuate  slowly 
with  height  until  a  distance  comparable  to  the  soliton  scale  is  achieved,  beyond  which  the  fields 
vanish  more  rapidly  .  Note  that  the  magnetic  fields  are  only  reduced  by  a  factor  of  2  at  500  m 
altitude.  This  could  be  contrasted  with  deep-ocean  internal  wave  magnetic  fields,  which  attenu¬ 
ate  exponentially  with  height  with  a  wavelength-dependent  e -folding  scale  that  does  not  exceed  a 
few  hundred  meters.  The  difference  in  behavior  is  due  to  the  localized,  dipolar  nature  of  the 
internal  soliton  sources,  while  the  deep-ocean  internal  waves  are  due  to  an  extended  and  random 
source. 

While  these  results  show  the  vector  behavior  of  internal  soliton  magnetic  fields,  most  air¬ 
borne  magnetometers  measure  only  the  scalar  magnitude  rather  than  the  actual  vector  com¬ 
ponents.  The  anomaly  seen  by  an  absolute  intensity  instrument  is  equivalent  to  the  projection  of 
the  induced  magnetic  field  onto  the  local  geomagnetic  induction  vector  because  the  latter  far 
exceeds  the  former  in  size  and  will  dominate  the  measurements.  The  total  field  anomaly  is  given 


where  B  is  the  internal  solilon  magnetic  field  (62)  (63)  and  (66)  (67)  and  F  is  the  geomagnetic 
field.  As  always,  there  are  two  limiting  cases  to  be  considered  a  soliion  traveling  to  the  north 
in  which  (62)“(63)  contribute,  and  soliton  moving  to  the  east,  in  which  only  induction  by  the 
vertical  force  is  important.  Figures  9  and  10  illustrate  these  two  situations.  Since  the  geomag¬ 
netic  field  is  nearlj  horizontal  and  north-directed  at  the  latitude  of  the  Andaman  Sea.  the  north- 
mo\ing  soliton  pr(>duces  a  sizable  anomaly  with  a  peak-to-peak  amplitude  of  3  .o  nT  at  the  sea 
surface.  By  wa\  of  contrast,  an  east-moving  soliton  gives  an  amplitude  of  only  about  0.2  nT 
The  remaining  behavior,  and  especially  the  attenuation  with  height,  is  like  that  in  Figures  7  and 
h 

Figures  11  and  12  illustrate  the  vector  magnetic  components  for  an  internal  soliton  in  the 
Gulf  of  California  using  the  hydrodynamic  parameters  of  Table  1  and  moving  to  the  north,  as  in 
Figures  7  and  The  peak-to-peak  amplitudes  are  about  1.5  iiT  at  the  sea  surface,  attenuating 
slowly  with  height  until  a  distance  comparable  to  the  soliion  scale  L  is  reached.  The  horizontal 
scale  of  the  magnetic  field  reflects  the  300  m  width  of  the  model  soliton.  Figures  13  and  14  show 
the  total  field  anomalies  obtained  from  (79)  for  solitons  traveling  to  the  north  and  east.  Al  this 
mid-latitude  site  the  horizontal  and  vertical  geomagnetic  fields  are  comparable  m  magnitude,  and 
only  shape  differences  are  seen  between  the  anomalies  for  the  two  directions  of  soliion  propaga¬ 
tion. 

Figures  15  and  16  show  the  horizontal  and  vertical  magnetic  components  for  the  Mas¬ 
sachusetts  Bay  case  in  Table  1.  again  for  a  north-moving  disturbance.  Due  to  the  very  small 
(60  m)  scale  of  the  soliton.  attenuation  for  heights  above  100  m  is  severe.  The  surface  peak-tc- 
peak  amplitude  is  about  0  5  nT  for  both  the  horizontal  and  vertical  magnetic  fields  Figures  17 
and  18  give  the  total  field  anomalies  for  solitons  propagating  to  the  north  and  east  respect  ivcK 
The  behavior  is  similar  to  that  of  the  vector  components  and  reflects  the  mid-latitude  equiparti- 
tion  of  the  geomagnetic  field  into  horizontal  and  vertical  parts. 
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Figures  19“21  give  the  corresponding  seafloor  niagneiic  fields  for  the  Andaman  Sea.  Gulf  of 
California  and  Massachusetts  Bay  models  using  (68)“  (71)  and  assuming  northward  soliton  propa¬ 
gation  In  each  case,  the  position  x  is  arbilrariK  set  tc-  zero  and  the  time  coordinate  has  been 
reversed  to  reflect  the  behavior  that  would  actually  be  observed  by  a  fixed  observer.  The  fields 
are  comparable  in  size  to  the  sea  surface  values,  while  the  time  scales  reflect  the  value  of  L  c  . 
where  L  is  the  soliton  scale  and  c  the  soliton  phase  velocity.  These  fields  are  quite  large,  espe- 
ciall\  when  compared  to  the  ambient  ionospheric  values  at  the  periods  of  interest,  and  could  be 
expected  to  contaminate  seafloor  geophysical  data  quite  appreciably. 

In  closing,  two  additional  observations  must  be  made.  First,  these  magnetic  anomalies  are 
nearly  static  in  behavior  and  induction  in  either  the  ocean  or  earth  is  negligible.  While  a  fixed 
seafloor  conductivity  of  0  OOo  S  m  was  assumed,  a  change  in  this  quantity  by  an  order  of  magni¬ 
tude  does  not  alter  the  magnetic  field  by  more  than  a  few  percent.  The  value  of  seawater  conduc¬ 
tivity  (fixed  at  3.2  S  m)  is  onh  important  in  its  effect  on  the  magnitude  of  the  source  currents 

(74)~(77)  with  (78).  and  is  well-known  in  any  case.  In  a  similar  vein,  the  frequency  in  (35)“ 

(42)  would  ordinarily  bo  replaced  by  the  frequency-of-encounter  kt  for  a  sensiir  moving  at  a  velo¬ 
city  f  which  is  much  larger  than  the  phase  velocity  c  of  the  soliton.  This  suggests  that  i  should 
replace  c  in  .4,  in  (62|~(63)  and  (66)“{71).  Since  induction  is  a  small  effect,  this  also  does  not 
alter  the  results  in  Figures  1  18  very  much.  This  suggests  that  the  Green  functions  (55)“(58) 
could  be  replaced  with  the  whole  space  values  without  introducing  a  substantial  amount  of  error, 
and  could  yield  analytic  values  for  the  magnetic  fields. 

Finally,  the  models  presented  in  this  report  are  for  an  isolated  soliton  of  depression.  In  the 

ocean,  internal  solitons  occur  in  packets  consisting  of  a  series  of  rank-ordered  solitons  whose 

wavelengths  and  amplitudes  decrease  from  its  leading  to  its  trailing  edge.  Fortunately,  the 
Maxwell  equations  are  linear,  and  a  real  internal  wave  packet  could  be  modeled  by  superimposing 
the  results  of  this  report  with  suitable  spatial  offsets  and  changes  in  the  soliton  parameters. 
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e  1  Internal  Soltlon  Hydrodxnamic 
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Figure  1.  Schematic  diagram  of  internal  wave  packet  generation  a?  a  lee  wave  during  ebh  tid 
(a),  followed  bj  nonlinear  steepening  of  the  lee  wa'e  depression  (b)  after  the  tide  turns-  Th 
change  in  tidal  flow  direction  allows  the  packet  to  propagate  shorwards  and  over  thf  baihvmetri 
feature  Nonlinear  and  dispersive  effects  then  allow  undulations  to  develop  behind  the  front  in  (' 
and  (d!  (taken  from  Haur>  el  al..  1979). 
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Figure  2.  Sch"Ria:i‘  diagram  of  an  inierr.a!  soliton  on  a  two-layer  fluid  at  finite  depth  The  soli- 
ton  If  a  wavf  <.f  depressic  r.  when  th'-  upper  layer  ihickneff  h  ,  if  smaller  than  ihe  lower  la>ei 
thickness  fi ;  The  dashed  lines  are  locii  of  constant  water  velocity  (isotachs).  while  the  arrows 
indicate  the  magnitude  and  direction  of  water  flow.  A  small  surface  soliton  accompanies  tht 
interna!  one  and  surfacf  rips  are  observed  above  the  descending  liml 
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Figure  3.  The  component  of  the  magnetir  field  given  b>  (62)  for  the  Andaman  Sea  soliton 
parameter*,  in  Table  1.  a  seawater  conductivity  of  S.2  S/m.  and  a  seafloor  conductivitv  of 
0  OOii  S  m  The  values  at  the  sea-surface  (r  =  0)  and  altitudes  of  100,  500.  and  1000  m  are 
shown.  See  le.xt  for  details 
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Figure  4.  The  B,‘  componeni  of  the  magnetic  field  given  b>  (63)  for  the  Andaman  Sea  soliion 
parameierf  in  Table  1.  a  seawater  conductivity  of  5.2  S/m,  and  a  seafloor  conductivity  of 
0.005  S  m.  The  values  of  the  sea-surface  (z  =  0)  and  altitudes  of  100,  500,  and  1000  m  are 
shown.  See  text  for  details 
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Figure  5.  The  B/  component  of  the  magnetic  field  given  by  (66)  for  the  Andaman  Sea  soliton 
parameters  in  Table  1.  a  seawater  conductivity  of  5.2  S/m  and  a  seafloor  conductivity  of 
0.005  S  m  The  values  at  the  sea-surface  (r  =  0)  and  altitudes  of  100,  500.  and  1000  m  are 
shown.  See  text  for  details 
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Figure  9.  The  lot&l  magnetic  field  anomaly  for  the  Andaman  Sea  soliton  propagating  to  the 
north  The  size  of  the  anomaly  is  controlled  almost  entirely  by  the  horizontal  geomagnetic  field 
due  to  the  low  geographic  latitude  at  this  site. 
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Figure  13.  The  vertical  magnetic  field  anomaly  produced  by  ihe  Gulf  of  California  internal  soli 
ton  propagating  to  the  north 
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Figure  14.  The  total  magnetic  field  anomaly  produced  by  the  Gulf  of  California  internal  solit 
propagating  to  the  ea«i 
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Figure  16.  The  vertjcal  magnetic  Field  produced  b>  the  Massachusetts  Bav  internal  soliton  pro¬ 
pagating  to  the  north 
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Figure  19.  The  homontal  (lolid)  and  vei.ical  (dashed)  magnetic  fields  at  the  seafloor  for  the 
Andaman  Sea  internal  soliton  propagating  to  the  north.  The  sensor  position  is  assumed  fixed, 
and  the  time-evoloution  of  the  magnetic  field  is  shown. 
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Figure  20.  The  horizontal  (aolid)  and  vertical  (daahed)  magnetic  field*  at  the  seafloor  for  the 
Gulf  of  California  internal  soliton  with  northward  propagation. 
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Figure  31.  The  horitontti  (lolid)  and  vertical  (duhed)  mainetic  fields  at  the  seafloor  for  the 
Musaehuietts  Bay  internal  soliton  with  eastward  propagation. 


